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ABSTRACT: We review recent experimental tests of the gravitational inverse-square law 
and the wide variety of theoretical considerations that suggest the law may break down in 
experimentally accessible regions. 
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1 INTRODUCTION 

1.1 Background 

Gravitation was the first of the four known fundamental interactions to be under- 
stood quantitatively and the first "grand unification" in physics. Isaac Newton's 
Theory of Universal Gravitation connected terrestrial phenomena (the "falling 
apple") with astronomical observations (the "falling Moon" and Kepler's Laws). 
This theory stood virtually unchallenged until Albert Einstein developed his rel- 
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ativistic theory of gravitation in 1917. Since then, General Relativity has passed 

all experimental tests and is today the standard model of gravitation. Yet some 
three centuries after Newton, gravitation remains one of the most puzzling topics 
in physics. Recently, a completely unexpected and fundamentahy new gravita- 
tional property was discovered using distant Type la supernovae: the apparent 
acceleration of the Hubble expansion (1, 2), which is as yet unexplained. Fur- 
thermore, gravitation is not included, and in fact not includable, in the imposing 
quantum field theory that constitutes the standard model of particle physics. 

There is a broad consensus that the two standard models are incompatible. 
The strong, weak, and electromagnetic interactions are explained as results of 
the quantum exchange of virtual bosons, whereas the gravitational interaction is 
explained as a classical consequence of matter and energy curving spacetime. Be- 
cause quantum field theories cannot describe gravitation and General Relativity 
predicts an infinite spacetime curvature at the center of a black hole, neither of 
these two standard models is likely to be truly fundamental. 

Connecting gravity with the rest of physics is clearly the central challenge of 
fundamental physics, and for the first time we have a candidate theory (string or 
M-theory) that may unify gravitation with particle physics. But the remaining 
theoretical problems have focused attention on possible new phenomena that 
could show up as deviations from the familiar inverse-square law (ISL) of gravity, 
generally at length scales less than a few millimeters, but sometimes also at 
astronomical or even cosmological distances. We review these speculations in 
Section 2. 

Although it is conventionally assumed that the ISL should be valid for sepa- 
rations from infinity to roug hly the Planck length {Rp = yfGhf^ = 1.6 x lO'^s 
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m), until a few years ago this assumption had only been precisely tested for sep- 
arations ranging from the scale of the solar system down to a few millimeters. 
The reasons for this are obvious: On the one hand, there are no independently 
known mass distributions on length scales larger than the solar system, and on 
the other hand, it is difficult to get enough matter in close enough proximity 
to obtain a background-free gravitational signal at length scales smaller than 1 
mm. This contrasts strongly with Coulomb's Law (and its electroweak general- 
ization), which has been tested for separations down to 10~^* m in e'^e~ leptonic 
interactions at high-energy colliders (3). Although Coulomb's Law has not been 
experimentally verified at large length scales (relative to laboratory dimensions), 
a null-type laboratory measurement looking for effects of the galactic electro- 
magnetic vector potential, A, rules out deviations due to a finite photon mass for 
length scales up to ~ 2 x 10^° m (4). 

l.Ll Parameterizations 

Historically, experimental tests of Coulomb's and Newton's inverse-square laws 
were used to set limits on violations that, for gravity, took the form 

nr) = G^. (1) 

Prom the perspective of Gauss's Law, the exponent 2 is a purely geometrical 
effect of three space dimensions, so this parameterization was not well-motivated 
theoretically. Instead, it is now customary to interpret tests of the ISL as setting 
limits on an additional Yukawa contribution to the familiar contribution, 
which in the gravitational case creates a potential 

y(r) = -G^^[l + «e-Al , (2) 
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where a is a dimensionless strength parameter and A is a length scale or range. 
The Yukawa contribution is the static limit of an interaction due to the exchange 
of virtual bosons of mass rub = h/ (Ac), where mi, is the boson mass; the Yukawa 
form is also useful in other contexts (see Section 2.2.1). 

Some investigators (see, e.g., (5)) have considered the possibility that a nonzero 
graviton mass could lead to a "pure Yukawa" gravitational potential V{r) = 
—Gmim2e~^^^/r, recognizing that this phenomenological form does not have a 
well-defined theoretical foundation (see ref. (6) for a different approach to the 
implications of a nonzero gravitational mass). Others have considered power-law 
modifications to the ISL (7): 



where a at is a dimensionless constant and ro corresponds to a new length scale 
associated with a non-Newtonian process. Terms with N = 2 and A'^ = 3 may be 
generated by the simultaneous exchange of two massless scalar and two massless 
pseudoscalar particles, respectively (8, 9, 10), while N = 5 may be generated by 
the simultaneous exchange of two massless axions (11) or a massless neutrino- 
antineutrino pair (12). 

In this review, we focus on the parameterization of Equation 2; any experiment 
that detects a violation of the ISL will indicate a strength, a, and a length scale, 
A, that characterizes the violation. Once a violation is detected, it will become 
necessary to determine the functional form of the violation. The parameterization 
of Equation 2 has strong implications for experimental tests of the ISL. Any one 
test of the law necessarily covers a limited range of length scales. Suppose, for 
example, one performs a Keplerian test, comparing the orbits of two planets 
orbiting a common sun. Clearly, the test is insensitive to values of A much less 
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than the orbit radius of the inner planet. It is also insensitive to ranges of A much 

larger than the orbit radius of the outer planet because both planets simply feel 

a renormalized Newton constant G = G{1 + a). Consequently, a great variety 

of experiments is needed to effectively explore a wide variety of length scales. 

This contrasts with limits on Yukawa interactions from "equivalence principle" 

tests, where a single experimental result for a composition-dependent acceleration 

difference typically provides a constraint on a for values of A ranging from the 

length scale of the attractor to infinity (see, e.g., (13)). 

1.2 Scope of This Review 

This review concentrates on experimental tests of the ISL at length scales of mil- 
limeters or less, and on the wide range of theoretical developments suggesting 
that new phenomena may occur in this regime. We also discuss speculations 
about possible ISL violations at much larger length scales that could have im- 
portant cosmological implications. A extensive review of experimental results at 
longer length scales (14) appeared in 1999; we update it in Section 4.5 below. A 
review of extra "gravitational" dimensions, with emphasis on collider signatures, 
has recently appeared in this journal (15). A recent review of tests of the ISL 
from microns to centimeters is Reference (16). Our review is focused on work 
done since 1995 and should be current as of January 2003. An earlier review (13) 
covered spin-dependent forces that we do not consider here. 
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2 THEORETICAL SPECULATIONS 



2.1 Unifying Gravity with Particle Physics: Two Hierarchy Problems 

The two greatest triumphs of twentieth-century physics are general relativity and 
quantum mechanics. However, we do not currently know how to link these two 
theories, or how to do calculations consistently in situations where both gravity 
and quantum effects are important, such as near the Big Bang and the cores 
of black holes. Clearly general relativity must be contained in a more funda- 
mental quantum theory that would allow sensible calculations even in extreme 
conditions. However, attempts to quantize general relativity have been plagued 
with difficulties. Although one can construct an effective quantum field theory 
of gravity and particle physics that is sufficiently accurate for many applications, 
the theory is infamously "nonrenormalizable" or nonpredictive — an infinite num- 
ber of free parameters are needed to describe quantum effects at arbitrarily short 
distances to arbitrary precision. 

All known nongravitational physics is includable within the standard model of 
particle physics — a quantum field theory in which the weak and electromagnetic 
interactions are unified into a single framework known as the electroweak theory. 
Symmetry between the weak and electromagnetic interactions is manifest above a 
scale of roughly 100 GeV. This unification scale, where the electroweak symmetry 
is spontaneously broken, is known as the electroweak scale. The electroweak scale 
is set by a condensate of a scalar field known as the Higgs field that has a negative 
mass-squared term of order (100 GeV)^ in its potential. All three forces of the 
standard model, the electromagnetic, weak, and strong interactions, are similarly 
unifiable into a simple group with a single coupling at the fantastically high energy 
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scale of 10^^ GeV. This "grand" unified theory (GUT) explains the quantization 

of electric charge and, provided there exists a new symmetry between fermions 

and bosons known as supersymmetry, predicts the observed value for the relative 

strengths of the weak and electromagnetic couplings. But supersymmetry has 

not yet been observed in nature and, if present, must be spontaneously broken. 

Supersymmetry and GUTs are reviewed in References (17, 18, 19, 20, 21, 22, 23, 

24, 25, 26, 27, 28). 

Intriguingly, the Planck scale, Mp = \/Uc[G^ at which quantum-gravity effects 
must become important, Mpc? = 1.2 x 10^^ GeV, is rather close to the apparent 
unification scale of the other forces. This hints that all belong together in a 
unified framework containing a fundamental scale of order Mp. Motivated by 
GUTs, the conventional view is that the phenomenal weakness of gravity at 
accessible energies — 10^^ times weaker than the other forces at the electroweak 
scale — is due to the small masses of observed particles relative to Mp. 

In the standard model, particle masses derive from the Higgs condensate. The 
tremendous discrepancy between the scale of this condensate and the presumed 
fundamental scale of physics is known as the gauge- hierarchy problem. In the 
minimal standard model, the smallness of the Higgs mass-squared parameter rel- 
ative to the GUT or Planck scales violates a principle known as "naturalness" — 
renormalized values of parameters that receive large quantum corrections should 
not be much smaller than the size of the corrections. The Higgs mass squared 
receives corrections proportional to the cutoff or maximum scale of validity of 
the theory. Naturalness would therefore demand that to describe physics at en- 
ergies higher than about 1 TeV, the standard model should be contained within 
a more fundamental theory in which the quantum corrections to the Higgs mass 
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are suppressed. An example of such a theory is a super symmetric extension 

of the standard model. In theories with spontaneously or softly broken super- 
symmetry, the quantum corrections to scalar masses are proportional to the 
supersymmetry-breaking scale. Provided the supersymmetry-breaking scale is 
of order 100 GeV, the electroweak scale is natural, and the hierarchy question is 
why the supersymmetry-breaking scale is so small compared with Mp. This latter 
problem is theoretically tractable; in many super symmetric models, the scale of 
supersymmetry breaking is proportional to exponentially small, nonperturbative 
quantum effects (29, 30). 

A second, and much bigger, hierarchy problem is known as the cosmological- 
constant problem. The strong observational evidence (1, 2) that the expansion 
of the universe is accelerating can be explained by a nonvanishing cosmological 
constant. The concordance of cosmological data indicates (31) that the universe is 
filled with a vacuum-energy density 

Pvac ~ 0-7/3c) where is the critical density 
3i?^c^/(87rG) and H is the present value of the Hubble constant. This gives 
Pvac ~ 4 keV/cm^, which corresponds to an energy scale ^ {TicYpvac ~ 2 meV or 
a length scale ^ {he)/ p 

vac ~ 100 yum. Such a small energy density is particularly 
puzzling because the quantum corrections to the vacuum energy density from 
particle physics scale as the fourth power of the cutoff of the effective theory. 
Such a cutoff might be provided by new physics in the gravitational sector. The 
energy scale of new gravitational physics has been presumed to be around Mp, 
which would imply a cosmological constant 10^^° times larger than observed. 
The success of the particle physics standard model at collider energy scales is 
inconsistent with a cutoff lower than 1 TeV. Even a relatively low TeV cutoff gives 
a theoretical contribution to the cosmological constant that is some 10^^ times 
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larger than experiment. References (32) and (33) conjecture that this monstrous 

discrepancy could be eliminated with a much lower cutoff for the gravitational 

sector of the effective theory, around 1 meV, corresponding to new gravitational 

physics at a distance of about 100 /xm. The theoretical framework for such a low 

gravity scale is necessarily very speculative. However, just as the gauge hierarchy 

compels experimental exploration of the TeV scale, the cosmological-constant 

problem strongly motivates submillimeter-scale tests of gravity. 

General relativity itself gives indications that the theory of quantum gravity 
is radically different from a conventional quantum field theory. For instance, 
in theories of gravity, the concept of entropy must be generalized because en- 
tropy cannot be an extensive quantity scaling like volume. In fact, strong ev- 
idence favors an upper bound on the entropy of any region that scales as the 
surface area of the boundary of the region (34, 35, 36). A further conjecture, 
the "holographic principle," suggests that this entropy bound indicates that the 
fundamental degrees of freedom of a gravitational theory can actually be formu- 
lated in a lower-dimensional theory. Reference (37) reviews these ideas and their 
subsequent development. 

M-theory is a popular candidate for a theory of quantum gravity. This theory 
was called string theory when it was believed that its fundamental degrees of 
freedom were one-dimensional objects propagating in a 10-dimensional spacetime. 
Six of these dimensions were assumed to be rolled up into a compact manifold 
of size ~ Rp and unobservable. We now know that "string" theory necessarily 
contains many types of objects, known as "branes" or "p-branes," where p, the 
number of spatial dimensions of the p-brane, can be anywhere from to 9. This 
realization has revolutionized our understanding of string theory. Furthermore, 
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string theory is "dual," or physically equivalent as a quantum theory, to an 11- 

dimensional theory known as M-theory. There is much theoretical evidence that 

all known consistent string theories, as well as 11-dimensional supergravity, are 

just weakly coupled limits in different vacua of a single theory of quantum gravity. 

Extra dimensions might seem to contradict the holographic assertion that the 
fundamental theory is actually Zower-dimensional. However, as comprehensively 
reviewed in Reference (38), the discovery that string theory on certain spacetimes 
with n noncompact dimensions is dual to a nongravitational gauge theory with 
n — 1 dimensions provides additional theoretical evidence for holography, as well 
as for string theory. Strings, M-theory, p-branes, and duality have been reviewed 
extensively (39, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55) and are the 
subject of several excellent textbooks (40, 41). 

Until recently, it was believed that experimental verification of a theory of 
quantum gravity was out of the question, due to the impossibly short distance 
scale at which quantum gravitational effects are known to be important. Fur- 
thermore, string theory contains a stupendous number of vacua — with no known 
principle for selecting the one we should live in — and so appears to have limited 
predictive power. Its chief phenomenological success to date is that in many of 
these vacua, the low-energy effective theory approximately resembles our world, 
containing the fields of the standard model and gravity propagating in four large 
dimensions. A major unsolved difficulty is that all known vacua are super sym- 
metric, although there are a variety of conceivable ways for the supersymmetry 
to be broken by a small amount. 

As we discuss below, although string theory makes no unique prediction, all 
known ways of rendering our observations compatible with string theory lead 
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to new, dramatic signals in feasible experiments. In particular, the discovery of 

branes has led to new possibilities for explaining the gauge hierarchy and the 

cosmological constant. Many of these can be tested in measurements of gravity 

at submillimeter scales, or in searches for small deviations from general relativity 

at longer distances. 

2.2 Extra Dimensions and TeV-scale Unification of Gravity 
2.2.1 "Large" extra dimensions 

It is usually assumed that the Planck scale is an actual physical scale, as is 
the weak scale, and that the gauge-hierarchy problem is to explain the origin 
of two vastly disparate scales. However, Arkani-Hamed, Dimopoulos and Dvali 
(ADD) (56) have proposed an alternative explanation for the weakness of grav- 
ity that has stimulated much theoretical and experimental work (see reviews in 
(15, 86, 87, 88, 89, 90)). Arkani-Hamed et al. conjecture that gravity is weak, 
not because the fundamental scale is high but because gravity can propagate in 
new dimensions less than a millimeter in size. Such "large" new dimensions are 
not seen by the standard-model particles because these are confined to a three- 
dimensional subspace of the higher-dimensional theory. Such a framework can 
be accommodated in string theory (57) . A type of p-brane known as a Dp-brane 
does have gauge and other degrees of freedom as light excitations that are con- 
fined to the brane. If the standard-model particles are all confined to such a 
D3-brane, we will not sense additional dimensions except via their modification 
of the gravitational force law. 

The hierarchy problem can be reformulated in this framework. One can assume 
that the fundamental scale M* is of order 1 TeV (58). There is then no hierarchy 
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between the weak scale and and no gauge-hierarchy problem. If there are n 

new dimensions, the higher-dimensional Newton's constant G(4+„) can be taken 

to be 

where 5'(2+„) is the area of a unit (2 + n)-sphere, 

At sufficiently short distances, the gravitational force at a separation r would 
be proportional to G(4+„)/r^"'"". To reconcile this with the 1/r^ force law ob- 
served at long distances, Arkani-Hamed et al. take the n new dimensions to be 
compact. At distances that are long compared with the compactification scale, 
the gravitational flux spreads out evenly over the new dimensions and is greatly 
diluted. Using Gauss's Law, one finds that for n new dimensions with radius i?*, 
compactified on a torus, the effective Newton's constant at long distances is 



Tic 



n 



1 n 



t'TJ, 



_M*ci 

Here is the volume of the n-torus, (27ri?*)". The relationship between R^ and 
M^, for other geometries may be found simply by using the appropriate formula 
for the volume. 

The hierarchy problem is then transmuted into the problem of explaining the 
size of the new dimensions, which are much larger than the fundamental scale. 
There are several proposals for stable compactifications of new dimensions that 
are naturally exponentially large (59, 60, 61, 62, 63). 

To test the ADD proposal directly, one should probe the ISL at a distance scale 
on the order of i?*. Compact new dimensions will appear as new Yukawa- type 
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forces, of range R^, produced by the exchange of massive spin-2 particles caUed 
Kaluza-Klein (KK) gravitons (64, 65, 66). To see this, note that the components 
of the graviton momenta in the compact dimensions must be quantized. For 
instance, compactification of a flat fifth dimension on a circle of radius R would 
impose the condition on P5, the fifth component of the graviton momentum, 
P5 = j^/Rj where j is an integer. The dispersion relation for a massless particle 
in five Lorentz-invariant dimensions is 

E' = Y.c'Pi + c'Pi. (7) 

i=0 

Comparing this with the four-dimensional massive dispersion relation 

E^ = j2c'Pi + c''M\ (8) 

1=0 

we see that the fifth component of the momentum appears as a four-dimensional 
mass term. A five-dimensional graviton thus appears as an infinite number of 
new massive spin-2 particles. For a fiat new dimension compactified on a circle 
of radius R, the mass nij of the jth KK mode is nij = j%/{Rc) with j = 1, 2, .... 

In factorizable geometries (whose spacetimes are simply products of a four- 
dimensional spacetime with an independent n-dimensional compact space), the 
squared wave functions of the KK modes are uniform in the new dimensions. Low- 
energy effective-field theory analyses of the KK modes and their couplings (67, 
68, 69, 70) show that higher-dimensional general coordinate invariance constrains 
this effective theory. Even at distances less than R, KK mode exchange will not 
violate the equivalence principle. The leading terms in an expansion in 1/M* 
contain a universal coupling of each graviton KK mode G'^^ to the stress tensor 
of form 




(9) 
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that is, each KK mode simply couples to the stress tensor in the same manner 

as the graviton. To compute the correction to the ISL for nonrelativistic sources 

at long distances, it suffices to consider the correction to the potential from the 

exchange of the lightest KK gravitons. The propagators for the KK states may 

be found in References (67, 68, 69, 70). 

For n new dimensions compactified on a flat torus, with the same radius i?* 
for each dimension, the lowest-lying KK mode has multiplicity 2n and Compton 
wavelength R^. Direct searches for such new dimensions would observe such KK 
gravitons via the contribution of their lowest-lying modes to the Yukawa potential 
of Equation 2, giving a = 8n/3 and A = ii*. A factor of 4/3 occurs in a because 
a massive spin-2 particle has five polarization states, and the longitudinal mode 
does not decouple from a nonrelativistic source ^ . Other compact geometries will 
give similar effects, although the value of a is quite model-dependent. 

Assuming all new dimensions are compactified on a torus of radius i?*, and 
= 1 TeV, Equation 6 gives 

R* ilQ-^^+f cm . 

TT 

The case n = 1, i?* = 3 x 10^^ m, is clearly ruled out. The case n = 2, R^, = 0.3 

mm, is inconsistent with the results of Reference (71). This case is even more 

strongly constrained by the observation of the neutrinos from supernova 1987A 

(72, 73, 74, 75, 76). Gravitational radiation into the extra dimensions would 

rapidly cool the supernova before the neutrinos could be emitted, imposing a con- 
^Note that References (81) and (82) included a contribution from a massless "radion" (gravi- 
tational scalar) in their Newtonian potential, and the radion KK modes in the Yukawa potential, 
leading to a different value for a. We discuss the radion and why it should be massive later in 
this section. 
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straint i?* < 0.7 //m. The extra gravitational degrees of freedom also necessarily 

spoil the successful calculations of big-bang nucleosynthesis unless < 2 /xm, 

and the decay of the KK modes would add a diffuse background of cosmological 

gamma rays whose non-observation implies i?* < 0.05 /xm (77). For n > 3, i?* is 

less than about a nanometer, which is still allowed by astrophysics, cosmology, 

and direct searches. 

It might, therefore, seem that direct observation of the new dimensions in 
ISL tests is out of the question. This conclusion is false. Astrophysical and 
cosmological bounds are still consistent with a single extra dimension of size 1 
mm — in such a scenario the hierarchy problem might be solved via the existence 
of several more much smaller new dimensions (78) . Furthermore, as discussed in 
the next section, it is easy to alter Equation 6 and the predictions for higher- 
dimensional graviton emission. Finally, there is a strong argument that the ADD 
proposal should modify the ISL at a scale of order hMp/{cM^). 

In theories of gravity, the geometry of spacetime is dynamical and can fluctuate. 
In particular, the radius of new dimensions can fluctuate independently at each 
point in our four-dimensional spacetime. Thus, low-energy effective theories of 
compact extra dimensions inevitably contain spin-0 fields parameterizing the radii 
of the new dimensions. If the size of the new dimensions is not determined 
by dynamics, then the linear combination of these fields that determines the 
extra dimensional volume is a massless Brans-Dicke scalar (79) with gravitational 
strength coupling, known as the "radion." A massless radion is decisively ruled 
out by tests of general relativity (80). Stabilization of the volume of the extra 
dimensions is equivalent to a massive radion. Since, with a low fundamental scale, 
the effective potential for the radion should not be much larger than 0{Mf), 
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and its couplings are proportional to Gn, the radion mass squared should be 
lighter than 0{GnM^). The radion will mediate a new, gravitational strength 
force, with a = n/(n + 2) ((85); G. Giudice, R. Rattazzi, N. Kaloper, private 
communications). In many cases, the radion is the lightest state associated with 
new dimensions. For M* less than a few TeV, its range should be longer than 
of order 100 jim.. Even for relatively "small" new dimensions, with size of order 
an inverse TeV, the radion will, under certain assumptions, have a Compton 
wavelength in the vicinity of 100 jiia (83, 84). 

2.2.2 Warped extra dimensions 

The previous discussion assumed the metric for the new dimensions is factor- 
izable. However, the most general metric exhibiting four-dimensional Poincare 

invariance is a "warped product," 



where the are the coordinates of the new dimensions, and / and g are general 
functions of those coordinates. Solving the higher-dimensional Einstein equations 
for a spacetime with an embedded brane with nonvanishing tension typically 
requires warping. The "warp factor" /(^i) may be thought of as a ^-dependent 
gravitational redshift factor that leads to a potential term in the graviton wave 
equation. This potential can have a dramatic effect on the ^ dependence of the 
wave functions of the graviton, the graviton KK modes, and the radion. 

Randall &: Sundrum (91) (RS-I) noted that a large hierarchy can be obtained 
with a single small new dimension if the metric takes the form 



(10) 



(11) 
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where ^ is a coordinate living on the interval [0, vr], A: is a constant, and Tc is 
the compactification scale. This is just the metric for a slice of five-dimensional 
anti-deSitter space (maximally symmetric spacetime with constant negative cur- 
vature). It is also a solution to the five-dimensional Einstein equations with 
five-dimensional Newton's constant 1/M^ if there is a negative cosmological con- 
stant of size A = — 24M^fc^, and if 3-branes are located at ^ = and ^ = tt with 
tensions ±24M^k. A negative-tension brane seems unphysical, but such bizarre 
objects can be constructed in string theory, provided the spaces on each side of 
the brane are identified with each other, that is, the brane represents a bound- 
ary condition on the edge of space. For large kvc, most of the extradimensional 
volume of this space is near the positive-tension brane at ^ = 0. 

To study the long-distance behavior of gravity in such a spacetime, one exam- 
ines the behavior of small fluctuations of this metric of the form 



Here h^i, is the four-dimcnsional graviton. Plugging this metric into Einstein's 
equations and linearizing in h, one finds h is a zero mode, or massless solution to 
the equations of motion, whose wave function in the compact dimension simply 
follows the warp factor e~^''^''^. Thus, there is a massless four-dimensional gravi- 
ton that is localized about the brane at ^ = and exponentially weakly coupled 
to matter on the brane at ^ = tt. If we further hypothesize that the latter brane 
is where the standard model lives, the weakness of gravity is explained for a mod- 
erate value of krc ~ 12. Both k and can be of the same order of magnitude as 
the fundamental scale, and so there is no large hierarchy in the input parameters. 

As in the ADD case, the RS-I model has a radion parameterizing the com- 
pactification scale. Goldberger &; Wise (92) have shown that krc in the desired 



(12) 
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range can naturally be stabilized without large dimensionless inputs if the theory 

contains a massive scalar that lives in the bulk and has source terms localized on 
the branes. The radion then acquires a large mass of order 100 GeV. The curva- 
ture in the extra dimension has a huge effect on the KK graviton spectrum and 
couplings. The lightest KK modes have masses in the TeV region and large wave 
functions near our brane, and therefore 0{1) couplings to ordinary matter. This 
model has unusual experimental signatures at colliders (15) but is not testable 
with feasible probes of the ISL. 

The RS-I model teaches us that warping can have significant effects on the 
phenomenology of the new dimensions. The coupling strength and masses of 
both the KK modes and the radion can be altered, and the graviton can be 
localized, or bound to a brane. Furthermore, warping is a generic phenomenon 
that should also occur in the ADD scenario. Even a very small amount of warping 
can greatly alter the coupling of the zero-mode graviton to our brane, which 
makes this coupling either much stronger or much weaker than for the case of 
flat extra dimensions (93), altering the relation of Equation 6. Even in the case 
of M^, = 1 TeV and n = 2, with a very small amount of warping, the masses of 
the lightest KK modes can be either higher or lower than the inverse-millimeter 
scale predicted by the unwarped case. 

2. 3 Infinite- Volume Extra Dimensions 

In a second paper (94) , Randall &: Sundrum (RS-II) explored the phenomenology 
of a graviton zero mode that is localized about a 3-brane embedded in a noncom- 
pact, infinite extra dimension. They found that five-dimensional gravity persists 
at all distance scales, with no gap in the KK spectrum, but at long distances the 
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1/r^ force, mediated by the zero mode bound to the brane, dominates, and the 

extra dimension can be unobservable at low energy. A simple model of this effect 

is given by the metric 

ds^ = e-'^^\\^^dxi'dx'' + dz^, (13) 

where z, the coordinate of the fifth dimension, is noncompact. This metric, which 
represents two slices of anti-deSitter space glued together at z = 0, also solves 
Einstein's equations, given a negative bulk cosmological constant — 24M^A;^, and 
a single 3-brane at 2; = of positive tension 24M^fc. The total gravitational 
potential between two masses mi and m2 separated by a distance r on the brane 
may be found by summing up the contributions of the bound-state mode and the 
continuum KK spectrum, which, for distance scales longer than l/fc, gives 

V(r)=G.r^(l^^) (14) 

with Gn = h'^k/M^. The experimental upper bound on 1/k from N = 3 terms in 
Equation 3 has not been explicitly computed but should be similar to the bound 
on the radius of an extra dimension. Therefore M* must be larger than about 
10^ GeV, and there is still a gauge hierarchy. With two or more infinite new 
dimensions, and a graviton confined to our 3-brane, it is possible to lower 
to 1 TeV (95). In such a scenario, the weakness of gravity is due to the zero- 
mode graviton wave function spreading over the extra dimensions, as in the ADD 
proposal, but the width of the wave function is set by the curvature scale rather 
than by the size of the dimension. Empirically, the main distinction between 
such weak localization and a large new dimension is that there is no gap in the 
KK spectrum and the ISL is modified by additional power-law corrections rather 
than by new Yukawa forces. 
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The RS-I explanation of the weakness of gravity — we hve on a brane, the 

graviton is confined to a different, parallel brane and its wave function here is 
small — can also be realized in infinite extra dimensions (95, 96). Lykken 8z Ran- 
dall studied such a configuration with a single extra dimension and concluded 
that the weakness of gravity could be explained without input of any large di- 
mensionless numbers. The chief test of their scenario would be strong emission 
of graviton KK modes at a TeV collider. The continuum of KK modes would 
modify the ISL, but their effect would only be significant for distances smaller 
than ~ 10 fm. 

2.4 Exchange Forces from Conjectured New Bosons 

Even if new dimensions are absent or small, the ISL can be modified at accessible 
distance scales by the exchange of new spin-0 or spin-1 bosons; spin-0 bosons 
would mediate an attractive Yukawa force while spin-1 bosons give a repulsive 
modification. Here we review some general considerations that apply to new 
bosons, and motivations for considering their existence. 

2.4.1 Scalars: general theoretical considerations 

In order for a scalar particle, (f), to exert a coherent force on matter, it must 
have a Yukawa coupling to electrons, to u, d, or s quarks, to the square of the 
gluon field strength, or to higher-dimension operators such as certain four-quark 
operators. The candidates of lowest dimension are 

y0ee, ^^d , ^^Puu , j^G^G'^'^'^ . (15) 

When embedded in the standard model, these all arise from dimension-5 op- 
erators, hence the common factor of 1//, where / has dimensions of mass. We 
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have assumed that aU chiral-symmetry-breaking operators should be proportional 

to fermion masses. With this assumption, and with all of the above operators 

present, the gluon coupling will dominate the scalar coupling to matter. Because 

the matrix element of in a nucleon is roughly the nucleon mass. Mat, such 

an interaction would lead to a Yukawa potential of the form given in Equation 2 

with A = ^/(m<^c), where is the scalar mass and a ~ M|,/(47r/^). 

An interaction {(f)/ f)G'^ produces radiative corrections to m^. In the standard 

model with cutoff A, one finds 

Sm^ ~ 4^ ^ "^"^ ■ ^^^^ 

The inequality expresses the requirement of naturalness. For / = Mp and = 
2 X 10~^ eV, corresponding to a Compton wavelength of 1 mm, naturalness 
implies A 5 TeV. This scale A approximately coincides with the scale at which 
naturalness of the electroweak-breaking sector demands new physics. A scalar 
coupled more weakly would correspond to a higher value for A. 

2.4.2 Forces from axion exchange 

A major loophole in the above arguments is that the interactions between matter 
and a new scalar may not arise from any of the operators in Equation 15, but 
rather from nonperturbative QCD effects. This is the case for the pseudoscalar 
axion invented to explain why strong interactions conserve CP to high precision. 
A pseudoscalar particle would normally not produce a Yukawa force between 
unpolarized bodies, but instantons in the presence of CP violation induce a 
scalar Yukawa coupling of the axion to matter that melts away above Aqcd- 
The softness of that coupling makes the radiative correction to the axion mass 
insignificant. However, a CP- violating scalar axion Yukawa coupling to matter 
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scales roughly as m„6QCD//a - &CiCD{muma)/{m„fT^), where < 5 MeV is 

the u quark mass, and 0qcd ^ 10~^(97) is the strong CP-violating angle. 

Thus, for an axion mass rua = 10~^ eV, the scalar axion coupling is at most 

about 10~^ times gravitational strength. ISL tests with unpolarized bodies probe 

the square of this coupling, so they are quite insensitive to the axion. On the 

other hand, monopole-dipole tests (98), which search for a CP- violating force 

between unpolarized and polarized bodies, are linear in the coupling and should 

be a more sensitive axion probe. 

2.4.3 Scalars: cosmological considerations 

A light, weakly interacting particle cannot decay or annihilate within a Hubble 
time, so its relic energy abundance must be equal to or less than that of the 
observed dark matter. However, the cosmology of scalars presents an important 
difficulty. A natural potential for a scalar in an effective theory below a cutoff A 
has the form V ~ A'^V{(p/f), where A ^ \Jvn^j ^ and V is an arbitrary function 
that is assumed to contain no large dimensionless numbers. If all scalar couplings 
are proportional to 1//, then the scalar lifetime is of order ATrf^/m^, essentially 
stable. If at a temperature T ~ A the thermal average of the scalar potential 
energy is {V) ~ T^, then the scalar field would have a large expectation value, 
(p ^ f. The infinite-wavelength component of this expectation value will be 
frozen until the Hubble scale is of order and will subsequently act like cold 

dark matter. Assuming the standard-model spectrum and standard cosmology 
for r < A (e.g., that the reheat temperature following inflation is above A), then 
an initial scalar energy density of at T = A implies a ratio today of the energy 
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in cold scalars to the energy in baryons of order 

^ ^ 2 X 10« f A) , (17) 

which is clearly unacceptable. 

Cosmology with light scalars can be made acceptable by invoking a very late 
stage of inflation with Hubble constant H less than or approximately m^. Then (p 
rapidly evolves to the minimum of its potential. Once inflation ends, the universe 
must reheat to a temperature Tr. However, the minimum of the scalar potential 
at Tr does not coincide with the minimum today, due to the tadpole generated 
by the interactions of Equation 15 at finite temperature. One must therefore 
check that coherent scalar oscillations are not regenerated during the reheating 
process after inflation. If reheating causes the minimum of the potential to change 
suddenly, relative to the oscillation time (of order 10~^^ s), then regeneration of 
the scalar condensate can be significant. We are almost completely ignorant of 
both the late inflationary mechanism and the timescale for reheating Ir, but a 
rough bound on may be estimated from the reheating temperature using the 
sudden inflaton-decay approximation 

,,.(./3)^-..,V3)(*)(/^. (18) 

For Tr > 10 MeV, which is necessary for standard big-bang nucleosynthesis, 
ti? ~ 3 X 10~^ s. Much higher reheat temperatures might be necessary to generate 
the baryon-number asymmetry. For example, a reheat temperature of about 100 
GeV corresponds to a reheat time of order ~ 3 x 10~^^ s. 

Provided this timescale is much longer than the scalar oscillation time fij {m^(?), 
the evolution of the minimum of the potential can take place adiabatically, in- 
jecting little energy into the coherent mode. The requirement of such a late stage 
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of inflation with acceptable reheating constrains theories of particle physics near 

the weak scale but does not rule out the existence of light scalar s. 



2.4.4 Bosons from hidden supersymmetric sectors 

As we discussed in Section 2, new physics is expected at the TeV scale. One 
candidate for this new physics is supersymmetry, which is expected in unified 
theories, and which can explain the gauge hierarchy. Unbroken supersymmetry 
predicts an unobserved degeneracy between fermions and bosons, so supersym- 
metry must be broken at a scale of 100 GeV or higher. The most popular scenario 
involves supersymmetry breaking at a scale of Ms ~ lO-*^^ GeV in a "hidden" sec- 
tor that couples to our visible world only via gravity and interactions of similar 
strength. The apparent scale of supersymmetry breaking in the visible world 
would then be of order Mg/Mp ~ 10^ GeV. In other scenarios, supersymmetry 
breaking is communicated to the visible world by the gauge forces of the stan- 
dard model, and the supersymmetry-breaking scale is as low as Ms ~ 10^ GeV. 
The supersymmetry-breaking scale is linked to m3/2, the mass of the gravitino 
(the spin- 1 supcrpartner of the graviton), through the relation m3/2 = M|/Mp. 
Well-motivated theoretical expectations for the gravitino mass range from 1 meV 
to 10^ GeV. In some scenarios (99, 100, 125, 126), the gravitino mass may be 
linked with the size of the cosmological constant inferred from the supernova 
observations and should be about 1 meV. 

If there are hidden sectors — particles coupled to the visible sector only via 
gravitational strength interactions — the apparent scale of supersymmetry break- 
ing in those sectors would typically be of order Scalar particles from those 
sectors could naturally have a mass in the meV range and mediate gravitational 
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strength forces with a range of about 100 //m. 

Note that the severe cosmological problems typical of light weakly coupled 
scalars discussed in the previous section do not necessarily occur for a scalar 
that is part of a hidden sector exhibiting super symmetry down to the meV scale. 
Such scalars might have a potential coming from 0{l) couplings to particles in 
the hidden sector, while maintaining a naturahy smah mass and gravitational- 
strength couplings to particles in the visible sector. These couplings will allow for 
the scalar field to relax to its minimum and for particle decay and annihilation. 

2.4.5 Forces from exchange of stringy bosons 

Supersymmetric hidden sectors are ubiquitous in string theory. All known ac- 
ceptable vacua of string theory are supersymmetric and contain a tremendous 
number of "moduli" — massless scalar fields whose expectation values set the pa- 
rameters of the effective theory. These moduli are extremely weakly coupled, 
with couplings inversely proportional to the fundamental scale. In order to give 
these fields a mass, it is necessary to break supersymmetry; however, moduli 
necessarily coTiplc weakly to the supersymmetry-breaking sector and, for a low 
supersymmetry-breaking scale, are expected to be extremely light. Current un- 
derstanding is inadequate to predict the moduli masses, but a rough estimate 
suggests these should be of order 777,3/2 (101). The best way to look for moduli is 
therefore to test the ISL at submillimeter distance scales. The couplings of the 
moduli in any given vacuum are computable, and so there are definite predic- 
tions. The best-understood scalar is the dilaton, a modulus that determines the 
strength of the gauge couplings. Its couplings to ordinary matter can be deter- 
mined and are nearly free of QCD uncertainties, so its discovery could provide a 
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genuine smoking gun for string theory (102). 



2.4.6 Forces from the exchange of weakly coupled vector bosons 

A new repulsive Yukawa interaction would be a signal for the exchange of a mas- 
sive spin-1 boson, presumably a gauge particle. In the ADD scenario, any gauge 
fields that propagate in the bulk of the new dimensions would have their couplings 
diluted by the same volume factor as the graviton and so would mediate a force 
with similar strength. Actually, since the gravitational force is also weakened by 
the smallness of the Mjv relative to M*, one would expect any such gauge forces 
to be stronger than gravity by a very large factor of (Af*/Afiv)^ ~ 10^-10^. This 
is acceptable if the range is substantially shorter than 1 mm (see Section 4.4.3). 
Gauge bosons could have a mass in an interesting range if the symmetry is broken 
via a scalar condensate on a brane. The resulting mass would be diluted by the 
bulk volume as well, and would naturally be in the range M*/ {VM^) ~ M^/Mp. 
For of order a few TeV, the range would be about 100 //m (72). If the sym- 
metry breaking occurs on the brane we live on, the gauge-boson couplings to 
standard-model matter could be substantially suppressed (105). 

Compactifications of string theory and other extradimensional theories often 
contain new massless spin-1 particles, known as graviphotons, that arise from 
components of the higher-dimensional graviton. These generally do not couple 
to ordinary light matter, but such bosons might acquire small masses and small, 
gravitational-strength couplings to ordinary matter, e.g., by mixing with other 
vector bosons (103, 104, 107, 108). Light spin-1 bosons do not suffer from the nat- 
uralness or cosmological difficulties of scalar particles, provided that they couple 
to conserved currents. However, spin-1 (and spin-0) boson exchange necessarily 
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"violates" the equivalence principle and the couplings of bosons with masses less 

than 1 neV are strongly constrained by the experiment of Reference (109). 



2.5 Attempts to Solve the C a smological- Constant Problem 

Comprehensive reviews of the cosmological-constant problem and the many at- 
tempts to solve it are available (110, 111, 112, 31, 113). Recent theoretical 
activity on this topic has been intense but is still inconclusive. Here we simply 
mention a few of the interesting recent proposals that imply modifications of the 
ISL at long distances. 

Beane (32) pointed out that in any local effective quantum field theory, natu- 
ralness would imply new gravitational physics at a distance scale of order 1 mm 
that would cut off shorter-distance contributions to the vacuum energy. Sun- 
drum (33) has speculated about the sort of effective theory that might do this. 
Sundrum proposed that the graviton is an extended object, with size of order 
1 mm, and has been exploring how to construct a natural and viable effective 
field theory from this picture (33, 114). It is still not clear how self-consistent 
this effective theory is, but it does have the great virtue of making a definite, 
testable experimental prediction — gravity should shut off below a distance scale 
of order 100 /xm. 

Many people have attempted to use extra dimensions to explain the small- 
ness of the cosmological constant, motivated by the alluring observation (115) 
that in higher-gravitational theories with branes, the four-dimensional vacuum 
energy or brane tension does not necessarily act as a source of four-dimensional 
gravity but can instead lead to curvature only in the new dimensions. So far 
no solved, consistent example actually yields a small cosmological constant in 
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the four-dimensional effective description without extreme fine tuning or other 

problematic features. 

Theories with branes and noncompact new dimensions allow another surpris- 
ing phenomenon known as quasilocalization of gravity (116, 117, 118, 119). In 
these theories, as in RS-II, long-distance gravity is higher-dimensional. However, 
there is no zero mode bound to our 3-brane. There is, instead, a metastable 
quasibound state that propagates four-dimensionally along the brane over times 
and distances that are short compared with some maximum scale. The ISL, 
and four-dimensional general relativity, will approximately apply from rmin to 
'"max, but not to arbitrarily long distances. The consistency of various theories 
of quasilocalization is still under debate and the theories themselves have been 
mutating rapidly. 

The holographic principle insinuates that a local description of a gravitational 
theory must break down somehow, because there are not enough degrees of free- 
dom to allow independent observables at different spacetime points. Several 
theorists have speculated that the breakdown of locality might even occur in a 
subtle way at astronomical or even longer distances, and that this might explain 
the size of the cosmological constant (120, 121, 122, 123, 124). In the Banks 
scenario (120, 123, 125, 126), supersymmetry ends up being broken at a scale of 
a few TeV by nonlocal effects due to the cosmological constant, leading to masses 
for the gravitino, dilaton, and other moduli of order 1 meV and deviations from 
the ISL at 100 /xm. 

Many of the above ideas share the possibility that there is some scale rmax 
beyond which Einstein gravity is modified. Modifying gravity at long distance 
allows a new approach to the cosmological constant. The observed acceleration 
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of the universe might be caused by a change in the behavior of gravity at the 

Hubble scale, instead of by dark energy (127). The fascinating prospect that 
the effective Newton's constant might be strongly scale-dependent at large dis- 
tance scales (gravity as a "high-pass spatial filter") leads to a new view of the 
cosmological-constant problem. Conventionally, it is assumed that the vacuum 
energy gravitates so weakly because, for some mysterious reason, this energy is 
actually very small. But if the strength of gravity depends on the wavelength of 
the source, it becomes credible that the vacuum energy is indeed very large but 
that it gravitates weakly because it is very smooth. Ideas along these lines have 
been pursued (128, 129, 130). 

References (131, 132, 133, 134, 135) present an intriguing assertion about the- 
ories of quasilocalization that may account for the acceleration of the universe. 
For any localized gravitational source, there exists a distance scale r* beyond 
which the graviton will acquire an extra polarization state that couples to the 
source so that the strength of gravity changes. This scale r* is a function of the 
gravitational radius of the source and Vmax', it decreases for less massive grav- 
itating objects. Dvali et al. (135) argue that ultraprecise measurements of the 
anomalous precession of the perihelion of planetary orbits can test models of 
quasilocalization that explain the cosmological acceleration. For instance, a 17- 
fold improvement in this measurement in the Earth-Moon system via lunar-laser 
ranging (LLR) would test a particular model in Reference (135). 
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3 EXPERIMENTAL CHALLENGES 



3. 1 Signals 

The dominant problem in testing gravitation at short length scales is the extreme 
weakness of gravity. This forces the experimenter to adopt designs that maximize 
the signal and minimize backgrounds and noise. For example, one could measure 
the force between spheres (136), between cylinders (137, 138), between a sphere 
and a plane (139, 140), or in planar geometry (71, 141). Clearly, at a given 
minimum separation, the signal from a short-range interaction, per unit test- 
body mass, is least for two spheres and greatest for two planes. 

The Yukawa force between two spheres of radii ri and r2 and masses mi and 
m2, whose centers are separated by s, is 

Fy = aG mim2$ f ^"j $ f ^"l f 1 + '—^ , (19) 



where ^{x) = 3(xcosha; — sinh a;)/a;^. For x S> 1, $(x) pa 3e^/{2x'^), whereas for 
X <C 1, <&(x) « 1. Therefore, for A <C r, the ratio of Yukawa to Newtonian forces 
for two spheres of radius r separated by a gap d is 

^«,^^fl + i^^,-^A. po) 



Fn 2r3 V 2r, 
The potential energy from a Yukawa interaction between a flat plate of area 
Ap, thickness tp, and density pp at distance d from an infinite plane of thickness 
t and density p is 

Vy = 2TraGpppX^Ap [l - e'^^/^] [l - e"*/^] e''^^^ , (21) 

if end effects are neglected. The corresponding force is 

Fy = 2TraGpppX^Ap \l - 6"*^ /^l [l - e"*/^] e""^/^ . (22) 
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In this case, for A much less than the thicknesses, the force ratio becomes 

^^a- e-'^A . (23) 
Fn tpt 

The potential energy of a Yukawa interaction between a sphere of radius r and 

mass m above an infinite plane of thickness t and density pp is 

Vy = iraGmpX^ $ (r / A) e"^/^ , (24) 

where s is the distance from the center of the sphere to the plane. The corre- 
sponding force is Fy = TraGmpX^{r / X)e~^/^ . In this case, for A <C r, the force 
ratio becomes 

^^«3^e-'^/\ (25) 
where d is the gap between the spherical surface and the plane. 

3.2 Noise Considerations 

Thermal noise in any oscillator sets a fundamental limit on the achievable sta- 
tistical error of its amplitude. A single-mode torsion oscillator subject to both 
velocity and internal damping obeys the equation 

T = ie + be + k{i + i4>)e , (26) 

where T is the applied torque, I the rotational inertia, 9 the angular deflection 
of the oscillator, and k the torsional spring constant of the suspension fiber. The 
velocity-damping coefficient h accounts for any losses due to viscous drag, eddy 
currents, etc., and the loss angle ^ accounts for internal friction of the suspension 
fiber. We compute the spectral density of thermal noise following Saulson's (150) 
treatment based on the fluctuation-dissipation theorem. The spectral density of 
torque noise power (per Hz) at frequency u is 

{Tl{u)) = 4kBT^{Z{u;)), (27) 
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where ks is Boltzmann's constant, T the absolute temperature, and Z = T/6 

the mechanical impedance. 

First consider the familiar case of pure velocity damping (6 > 0, = 0) where 

Z{lo) = iloj + b + k/ {iu). In this case, the spectral density of torque noise, 

{Tl{u;))=4kBT^-^ (28) 

{loq = is the free resonance frequency and Q = lujo/b the quality factor of 

the oscillator), is independent of frequency. The corresponding spectral density 
of angular-deflection noise in 9 is 

{OM) = Qj (a;2_ 0,2)2 + (^^^/Q)2 " (29) 

Note that the integral of Equation 29 over all / = uj/{2'k) is kBT/K, consistent 
with the equipartition theorem. The signal due to an external torque T is 

I^HI = J I , ] (30) 

so the signal-to-noise ratio in unit bandwidth has the form 
,H . ^frfL . , (31) 

where we have included a noise contribution (^^o) from the angular-deflection 
readout system. The signal is usually placed at a frequency u < to avoid 
attenuating the deflection amplitude 9 because of oscillator inertia. 

Now consider the case of pure internal damping (6 = 0, ^ > 0) where Z = 
iloj + K/{iLo) + K(f)/uj. In this case, the spectral density of thermal noise has a 
1// character, 

T. ,2 

{Tl{u))=AkBT-^, (32) 
where now Q = l/(p. The corresponding spectral density of thermal noise in the 



INVERSE-SQUARE LAW TESTS 35 
angular deflection is 

^•^^^^ ga;7(a;2-a;2)2 + (a;2/Q)2 • ^•^"'^ 
The signal-to-noise ratio in unit bandwidth is 

^ (34) 



so it is advantageous to boost the signal frequency above uiq until 6*^0 makes a 
significant contribution to the noise. 

S.S Backgrounds 

Electromagnetic interactions between the test bodies are the primary source of 
background signals and may easily dominate the feeble gravitational signal. In the 
following sections, we discuss the dominant electromagnetic background effects 
in ISL experiments. 

3.3.1 Electric potential differences and patch fields 

Electric charges residing on insulating or ungrounded test bodies are difficult to 
quantify, and Coulomb forces acting on such bodies can exceed their weights. For 
this reason, ISL tests typically employ conducting grounded test bodies. Even 
so, a variety of effects can give the test bodies different electric potentials. If 
dissimilar materials are used for the test bodies, a potential difference equal to 
the difference between the work functions of the two materials is present, typically 
of order 1 V. Even if the same material is used for both test bodies or the test 
bodies are both coated with the same material, such as gold, small differences 
in the contact potentials connecting the test bodies to ground can leave a net 
potential difference between the test bodies. With care, such contact potential 
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differences can be reduced to the level of a few mV (142)). 

Neglecting edge effects, the attractive electric force between a conducting plate 
with area A parallel to an infinite conducting plate is -Fe(c^) = eoAV^/{2(P), 
where d is the separation between the plates, V is the potential difference between 
the plates, and eo is the permitivity of free space. For 1-mm-thick plates with 
a density of 10 g/cm^, separated by 0.1 mm, Fe becomes as large as Fjv for 
a potential difference of 10 mV, and the electric force grows with decreasing 
separation whereas the Newtonian force is constant. 

Even if test bodies are at the same average potential, they experience a resid- 
ual electric interaction from patch fields — spatially varying microscopic electric 
potentials found on the surface of materials (143). Patch fields arise because 
different crystal planes of a given material have, in general, work functions (144) 
that can in extreme cases differ by as much as 1 V. To the extent that the surface 
is a mosaic of random microscopic crystal planes, local potential differences will 
occur with a scale size comparable to the size of the microcrystals. For example, 
the work functions of different planes of W crystals differ by 0.75 V. Gold is a 
good choice for test-body coating because the work functions of its crystal planes 
vary by only 0.16 V. Surface contaminants also contribute to the local variation 
of the electric potential, altering the local work function and providing sites for 
the trapping of electrical charge. In the limit that the patches are smaller than 
the separation, the patch field force (143) scales as 

3.3.2 Casimir Force 

Vacuum fluctuations of the electromagnetic field produce a fundamental back- 
ground to ISL tests at short length scales. The Casimir force (145) between 
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objects in close proximity may be viewed as arising either from the modification 

of the boundary conditions for zero-point electromagnetic modes or from the force 

between fluctuating atomic dipoles induced by the zero-point fields (146). The 

Casimir force can be quite large compared to the force of gravity. The Casimir 

force between two grounded, perfectly conducting, smooth, infinite planes at zero 

temperature, separated by a distance d, is attractive with a magnitude of 

A ~ 240^4 • ^^^^ 

For a 1-mm-thick plate of area A near an infinite plate of thickness 1 mm (again, 
both with density 10 g/cm^), Fc becomes equal to Fjv at a separation of d = 13 
Hm. 

Because precisely aligning two parallel planes is so difficult, experimenters usu- 
ally measure the force between a sphere (or spherical lens) and a plane. Assuming 
perfectly conducting, smooth bodies at zero temperature, the Casimir force is at- 
tractive with a magnitude of 

„ TT^Rhc , , 

where R is the radius of the sphere and d is the minimum separation between 
the surfaces of the sphere and plane. For a 1-mm-radius sphere near an infinite 
l-mm=thick plane (both with a density of 10 g/cm^), Fq becomes equal to F/v 
at a separation oi d = 2.5 /xm. 

The Casimir-force expressions in Equations 35 and 36 must be corrected for 
finite temperature, finite conductivity, and surface roughness (see below). All 
these corrections vary with the separation, d, making it difficult to distinguish a 
gravitational anomaly from an electrical effect. 
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3.3.3 Electrostatic shielding 

Fortunately, backgrounds from the Casimir force, electric potential differences, 
and patch-effect forces can be greatly reduced by using a moving attractor to 
modulate the signal on a stationary detector and placing a stationary, rigid, con- 
ducting membrane between the detector and the attractor. But this electrostatic 
shield places a practical lower limit of some tens of micrometers on the minimum 
attainable separation between the test bodies. 

3.3.4 Magnetic effects 

Microscopic particles of iron embedded in nominally nonmagnetic test bodies 
during their machining or handling, or in the bulk during smelting, can create 
local magnetic fields so small they are difficult to detect with standard magne- 
tometers, yet large enough to compete with gravitational forces. The magnetic 
force between two magnetically saturated iron particles 1 mm apart, each 10 //m 
in diameter, can be as large as 10~^ dynes, varying as the inverse fourth power of 
the distance between the particles. This is as large as the gravitational attraction 
between a 1-mm-thick Al plate with an area of 3 cm^ near an infinite Al plate 
that is 1 mm thick. Yet the magnetic field of such a particle is only 0.3 mGauss 
at a distance of 2 mm. 

Most ISL tests modulate the position of an attractor and detect the force this 
modulation produces on a detector. Even if the attractor has no ferromagnetic 
impurities, any magnetic field associated with the attractor modulation, e.g., 
from motor magnets or flowing currents, can couple to magnetic impurities in 
the detector. Experimenters typically measure the magnetic field associated with 
the modulation of the attractor and apply larger fields to find the response of the 
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detector. A variety of smaller magnetic background effects are associated with 

the magnetic susceptibilities of the test bodies. Standard magnetic shielding of 

the experimental apparatus is usually sufficient to reduce the ambient magnetic 

field to a level where the susceptibilities pose no problem. 

3.3.5 Other effects 

Modulation of the attractor position may introduce background effects that are 
not electromagnetic. The most obvious is a spurious mechanical coupling that 
transmits the motion of the attractor through the apparatus to the detector. 
These unwanted couplings can be reduced by multiple levels of vibration isolation 
and by experimental designs that force the signal frequency to differ from that 
of the attractor modulation. Experiments are performed in vacuum chambers to 
reduce coupling between the test bodies from background gas. 

3.4 Experimental Strategies 

ISL tests can be constructed as null experiments, partial-null experiments, or 
direct measurements. For example, Hoskins et al. (137) studied the force on a 
cylinder located inside a cylindrical shell. To the extent that the length-to-radius 
ratios of the cylinders are very large, this constitutes a null test because the 
Newtonian interaction between the cylinders gives no net force. Other null tests 
have used planar geometry; the Newtonian force between two parallel, infinite 
planes is independent of their separation. This basic idea, as discussed below, 
was exploited in Reference (141). An advantage of null experiments is that the 
apparatus does not need to handle signals with a wide dynamic range and the 
results are insensitive to instrumental nonlinearities and calibration uncertainties. 



40 ADELBERGER, HECKEL & NELSON 

Hoyle et al. (71) have reported a partial- null experiment in which the New- 
tonian signal was largely, but not completely, cancelled. As discussed below, 
the partial cancellation greatly reduced the required dynamic range of the in- 
strument, but Newtonian gravity still gave a very characteristic signal that was 
used to confirm that the instrument was performing properly. The form and 
magnitude of this signal provided constraints on new physics. 

Finally, Mitrofanov & Ponomareva (136) reported a direct experiment that 
compared the measured force beween two spheres as their separation was switched 
between two values. In this case, the results depended crucially on accurate 
measurement of the separations of the spheres and the forces between them. 

4 EXPERIMENTAL RESULTS 

4-1 Low- Frequency Torsion Oscillators 
4.1.1 The Washington experiment 

Hoyle et al. (71) of the University of Washington Eot-Wash group developed a 
"missing-mass" torsion balance (Figure 1), for testing the ISL at short ranges. 
The active component of the torsion pendulum was an aluminum ring with 10 
equally spaced holes bored into it. The pendulum was suspended above a copper 
attractor disk containing 10 similar holes. The attractor was rotated uniformly 
by a geared-down stepper motor. The test bodies in this instrument were the 
"missing" masses of the two sets of 10 holes. In the absence of the holes, the 
disk's gravity simply pulled directly down on the ring and did not exert a twist. 
But because of the holes, the ring experienced a torque that oscillated 10 times 
for every revolution of the disk — giving sinusoidal torques at lOw, 20a;, and 30a;, 
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where u was the attractor rotation frequency. This torque twisted the pendu- 
lum/suspension fiber and was measured by an autocollimator that reflected a 
laser beam twice from a plane mirror mounted on the pendulum. Placing the 
signals at high multiples of the disturbance frequency (the attractor rotation 
frequency) reduced many potential systematic errors. A tightly stretched 20-/xm- 
thick beryllium-copper electrostatic shield was interposed between the pendulum 
and the attractor to minimize electrostatic and molecular torques. The entire 
torsion pendulum, including the mirrors, was coated with gold and enclosed in 
a gold-coated housing to minimize electrostatic effects. The pendulum could not 
"see" the rotating attractor except for gravitational or magnetic couplings. Mag- 
netic couplings were minimized by machining the pendulum and attractor with 
nonmagnetic tools and by careful handling. 

The experiment was turned into a partial-null measurement by adding a sec- 
ond, thicker copper disk immediately below the upper attractor disk. This disk 
also had 10 holes, but they were rotated azimuthally with respect to the upper 
holes by 18° and their sizes were chosen to give a lOu torque that just cancelled 
the lOo; Newtonian torque from the upper attractor. On the other hand, a new 
short-range interaction would not be cancelled because the lower attractor disk 
was simply too far from the pendulum. The cancellation was exact for a sepa- 
ration (between the lower surface of the pendulum and the upper surface of the 
attractor) of about 2 mm. For smaller separations the contribution of the lower 
disk was too small to completely cancel the 10a; signal, and at larger separations 
the lower disk's contribution was too large (see Figure 2). 

Two slightly different instruments were used; both had 10-fold rotational sym- 
metry and differed mainly in the dimensions of the holes. In the first experiment. 
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the pendulum ring was 2.002 mm thick with 9.545-mm-diameter holes and a total 

hole "mass" of 3.972 g; in the second experiment, the ring thickness was 2.979 

mm with 6.375-mm-diameter holes having a total hole "mass" of 2.662 g. The 

resonant frequencies of the two pendulums, a;o/27r, were 2.50 mHz and 2.14 mHz, 

respectively; the fundamental Wu> signals were set at precisely ^ u>o and | u>o, 

respectively. In both cases, the 20a; and SOw harmonics were above the resonance. 

The observed spectral density of deflection noise was close to the thermal value 

given in Equation 33 for the observed Q factor of 1500 (see also Figure 3 below). 

4.1.2 Signal scaling relations 

The gravitational torque exerted on the pendulum by the rotating attractor is 
Tg{(f)) = —dV {(/))/ do, where V{(f)) is the gravitational potential energy of the at- 
tractor when the attractor is at angle (p, and 9 is the twist angle of the pendulum. 
For cylindrical holes, four of the six Newtonian torque integrals can be solved an- 
alytically but the remaining two must be evaluated numerically. Clearly, the 
Newtonian signal drops as the number of holes increases and their radii decrease 
because the long-range gravitational force tends to "average away" the holes. It 
also drops rapidly for separations much greater than the thickness of the upper 
attractor disk. Only three of the Yukawa torque integrals can be solved ana- 
lytically. However, when the Yukawa range. A, becomes much smaller than any 
of the relevant dimensions of the pendulum/attractor system, a simple scaling 
relation based on Equation 21 governs the signal and 

dA 

Ty oc aGppPaX'e-'/^ — , (37) 

where pp and pa are the densities of the pendulum and attractor, respectively; A 
is the Yukawa range; and A is the overlap area of the holes in the pendulum with 
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those of the attractor when the attractor angle is (p. 
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4.1.3 Backgrounds 

Hoyle et al. (71) found that the effects from spurious gravitational couplings, 
temperature fluctuations, variations in the tilt of the apparatus, and magnetic 
couplings were negligible compared with the statistical errors. Electrostatic cou- 
plings were negligible because the pendulum was almost completely enclosed by 
a gold-coated housing. The 20-)um-thick electrostatic shield was rigid to pre- 
vent secondary electrostatic couplings. The shield's lowest resonance was about 
1 kHz, and the attractor could produce a false electrostatic effect only by flexing 
the shield at a very high m = 10 mode. 

4.1.4 Alignment and calibration 

Although all submillimeter tests of the ISL face an alignment problem, it was 
especially important in this experiment because of the relatively large size of the 
pendulum (chosen to increase the sensitivity). Alignment was done in stages. 
First the pendulum ring was leveled by nulling its differential capacitance as the 
pendulum rotated above two plates installed in place of the electrostatic shield. 
The shield was then replaced, and the tilt of the entire apparatus was adjusted to 
minimize the pendulum-to-shield capacitance. To achieve horizontal alignment, 
the gravitational torque was measured as the horizontal position of the upper 
fiber suspension point was varied. Determining separations from mechanical or 
electrical contacts gave unreliable results, so the crucial separation between the 
pendulum and the electrostatic shield was determined from the electrical capac- 
itance. 
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The experimenters calculated the torque scale directly, using gravity. Two 

small aluminum spheres were placed in an opposing pair of the 10 holes of the 

torsion pendulum and two large bronze spheres, placed on an external turntable, 

were rotated uniformly around the instrument at a radius of 13.98 cm. Because 

this was close to the 16.76-cm radius (147) used in determining G and the ISL 

has been tested at this length scale (see Figure 4), the calibration torque could 

be computed to high accuracy. The torsion constant of the fiber was about 0.03 

dyne cm. 

4.1.5 Results 

Data were taken at pendulum/attractor separations down to 197 fj,m, where the 
minimum separation was limited by pendulum "bounce" from seismic distur- 
bances. The torque data, shown in Figure 2, were analyzed by fitting a potential 
of the form given in Equation 2 with a and A as free parameters and treat- 
ing the important experimental parameters (hole masses and dimensions, zero 
of the separation scale, torque calibration constant, etc.) as adjustable parame- 
ters constrained by their independently measured values. Hoyle et al. reported 
results from the first of the two experiments in Reference (71); the combined 
95%-confidence-level (CL) result of both experiments was given subsequently 
(148, 149) and is shown in Figure 5. 

The results exclude the scenario of two equal extra dimensions whose size gives 
a unification scale of M* = 1 TeV; this would imply an effective Yukawa inter- 
action with A = 0.3 mm and a = 16/3 if the extra dimensions are compactified 
as a torus. Because a > 16/3 is consistent with the data only for A < 130 fim, 
Equation 6 implies that > 1.7 TeV. A tighter bound on M* can be extracted 
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from the radion constraint, which, in the unwarped case where 1/3 < a < 3/4 

for 1 < n < 6, suggests that M* > 0(3 Tev). 

More interesting and general is the upper hmit placed on the size of the 
largest single extra dimension, assuming all other extra dimensions are signif- 
icantly smaller (71, 149). For toroidal compactification, this corresponds to the 
largest A consistent with a = 8/3, leading to an upper limit i?* < 155 /xm. Other 
compactification schemes necessarily give somewhat different limits. 

4-2 High- Frequency Torsion Oscillators 
4.2.1 The Colorado experiment 

The modern era of short-range ISL tests was initiated by Long et al. at the 
University of Colorado (151). Their apparatus, shown in Figure 6, used a planar 
null geometry. The attractor was a small 35 mm x 7 mm x 0.305 mm tungsten 
"diving board" that was driven vertically at 1 kHz in its second cantilever mode 
by a PZT (lead zirconate titanate) bimorph. The detector, situated below the 
diving board, was an unusual high-frequency compound torsion oscillator made 
from 0.195-mm-thick tungsten. It consisted of a double rectangle for which the 
fifth normal mode resonates at 1 kHz; in this mode, the smaller 11.455 mm x 
5.080 mm rectangle (the detector) and the larger rectangle (one end of which was 
connected to a detector mount) counter-rotated about the torsional axis, with 
the detector rectangle having the larger amplitude. The torsion oscillations were 
read out capacitively from the larger rectangle. The attractor was positioned so 
that its front end was aligned with the back edge of the detector rectangle and a 
long edge of the attractor was aligned above the detector torsion axis. A small 
electrostatic shield consisting of a 0.06-mm-thick sapphire plate coated with 100 
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nm of gold was suspended between the attractor and the detector. The attractor, 

detector, and electrostatic shield were mounted on separate vibration-isolation 

stacks to minimize any mechanical couplings and were aligned by displacing the 

elements and measuring the points of mechanical contact. 

In any null experiment, it is helpful to know the precise form of a signal of new 

physics. Long et al. slid away the electrostatic shield and applied a 1.5-V bias 

to the detector to give a large, attractive electrostatic force; this determined the 

phase of the signal that would be produced by a new, short-range interaction. 

4.2.2 Signal-to-noise considerations and calibration 

The spectral density of thermal-force noise in the multimode oscillator used in 
Reference (151) obeys a relation similar to Equation 32. The Colorado experi- 
menters operated on a resonance with a Q = 25, 000 so the readout noise was 
negligible. Data were taken with the attractor driven at the detector resonance 
as well as about 2 Hz below the resonance (see Figure 7). The mean values of 
the on-resonance and off-resonance data agreed within errors, but the standard 
deviation of the on-resonance data was about twice that of the off-resonance data. 
This is just what one would expect if the on-resonance data were dominated by 
thermal noise. Furthermore, the on-resonance signal did not change as the ge- 
ometry was varied. This ruled out the unlikely possibility that the observed null 
result came from a fortuitous cancellation of different effects, all of which should 
have different dependences on the geometry. The torsion oscillation scale was 
calibrated by assuming that the on-resonance signal was predominantly thermal. 
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4.2.3 Backgrounds 

Although a net signal was seen, it had the same magnitude on and off resonance 
and presumably was due to electronic pickup. No evidence was seen for an addi- 
tional, statistically significant background. Checks with exaggerated electostatic 
and magnetic effects showed that plausible electrostatic and magnetic couplings 
were well below the level of thermal noise. 

4.2.4 Results 

The null results from this experiment, taken at a separation of 108 //m, were 
turned into a(A) constraints using a maximum-likelihood technique. For various 
assumed values of A, the expected Yukawa force was calculated numerically 400 
times, each calculation using different values for experimental parameters that 
were allowed to vary within their measured ranges. A likelihood function con- 
structed from these calculations was used to extract 95%-CL limits on a(A). The 
results (141), shown in Figure 5, exclude a significant portion of the moduli forces 
predicted by Dimopoulos Sz Giudice (101). 

4-3 Microcantilevers 

4.3.1 The Stanford experiment 

Chiaverini et al. at Stanford (152, 153) recently reported a test of the ISL using 
the microcantilever apparatus shown schematically in Figure 8. This instrument 
was suited for the 10-/Ltm length scale but lacked the sensitivity to see gravity. The 
apparatus consisted of a silicon microcantilever with a 50 /Ltm x 50 /xm x 50 /xm 
gold test mass mounted on its free end. The cantilever had a spring constant of 
about 5 dyne/cm, and its displacement was read out with an optical-fiber inter- 
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ferometer. The microcantilever, which hung from a two-stage vibration-isolation 

system, oscillated vertically in its lowest flexural mode at a resonant frequency 
of ojQ 300 Hz . The microcantilever was mounted above an attractor consisting 
of five pairs of alternating 100 fiva xlOO /xm xl mm bars of gold and silicon. 
The attractor was oscillated horizontally underneath the cantilever at about 100 
Hz by a bimorph; the amplitude was chosen to effectively resonantly excite the 
cantilever at the third harmonic of the attractor drive frequency. The geometry 
was quite complicated; the third harmonic gravitational force on the cantilever 
depended sensitively and nonlinearly on the drive amplitude. An electrostatic 
shield consisting of a 3.0-//m-thick silicon nitride plate with 200 nm of gold evap- 
orated onto each side was placed between the cantilever and the attractor. Data 
were taken with the vertical separations between the cantilever and the attractor 
as small as 25 /xm. 

4.3.2 Signal-to-noise considerations 

The dominant noise source in the Stanford experiment was thermal noise in the 
cantilever, which was reduced by operating at about 10 K. The Q factors of the 
oscillating cantilevers in these measurements were typically about 1200. 

4.3.3 Calibration and alignment 

The cantilever spring constant k was found in two independent ways that agreed 
to within 10%: by assuming that when the cantilever was far from the attractor 
it was in thermal equilibrium with its surroundings, and by calculating k from 
the measured resonant frequency. The cantilever was aligned with respect to the 
attractor using magnetic forces. The cantilever's test mass had a thin nickel film 
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on one face, and the attractor was equipped with a zig-zag conducting path that 

followed the gold bars. When a current was run through the attractor, it placed 

a force on the cantilever that had half the frequency and phase of the expected 

gravitational signal but vastly greater amplitude. This force was used to align 

the apparatus. 

4.3.4 Backgrounds 

This experiment was limited by a spurious force about 10 times greater than 
the thermal detection limit. This force was clearly not fundamental, i.e., related 
to the mass distributions on the attractor, because the phase of signal did not 
behave as expected when the horizontal offset of the attractor oscillation was 
varied or as the attractor drive amplitude was changed. The most likely source 
of a spurious force is electrostatics; the cantilever was not metallized and so it 
could hold charge and the shield was observed to vibrate by a picometer or so. 
A potential on the cantilever of about 1 V would be sufficient to produce the 
observed force. Although thin nickel layers were incorporated into the test mass 
and attractor, the experimenters estimate that magnetic forces from the nickel 
(as well as from iron impurities in the gold) were too small to explain the observed 
background force. Vibrational coupling between the attractor and cantilever was 
minimized because the attractor was moved at right angles to the cantilever's 
flex. 

4.3.5 Results 

The experimenters saw a spurious (8.4ibl.4) x 10~^^ dyne force when the attractor 
and cantilever were at their closest separation of 25/xm. They assigned a 95%-CL 



50 ADELBERGER, HECKEL & NELSON 

upper limit on a Yukawa interaction by computing the minimum o; as a function 

of A that would correspond to this central value plus two standard deviations. 

Figure 9 shows their constraint, which rules out much of the parameter space 

expected from moduli exchange as computed in Reference (101). 

4-4 Casimir Force Experiments 

Early attempts to detect the Casimir force between metal surfaces (154) and 
dielectric surfaces (155, 156, 157, 158) had relatively large errors. Nonetheless, 
it was recognized (159, 160, 161) that such measurements provided the tightest 
constraints on new hypothetical particles with Compton wavelengths less than 
0.1 mm. In recent years, three groups have reported measurements of the Casimir 
force with relative errors of 1% to 5%. Although these experiments are orders 
of magnitude away from providing tests of the ISL, they do probe length scales 
from 20 nm to 10 /xm, where large effects may occur (see Section 2.4.4). 

4.4.1 Experimental methods 

The first of the recent experiments, performed by Lamoreaux at the University 
of Washington (139, 162), used a torsion balance to measure the force between a 
flat quartz plate and a spherical lens with a radius of 12.5 it 0.3 cm. Both surfaces 
were coated with 0.5 /im. of copper followed by 0.5 /xm of gold. A piezoelectric 
stack stepped the separation between the plate and lens from 12.3 fim to 0.6 /xm, 
at which point the servo system that held the torsion pendulum angle constant 
became unstable. The force scale was calibrated to 1% accuracy by measuring 
the servo response when a 300-mV potential difference was applied between the 
plate and lens at a large (ss 10 /xm) separation. The absolute separation between 
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the lens and plate was obtained by applying a potential difference between the 

two surfaces and fitting the measured force (for distances greater than 2 fim 
where the Casimir force was small) to the expected 1/d dependence, where d is 
the distance between the plate and lens. After subtracting the 1/d component 
from the force scans, the residual signals were fitted to the expected form for a 
Casimir force, and they agreed to within 5% (139, 162). 

Mohideen and collaborators at the University of California at Riverside re- 
ported a series of experiments that used an atomic-force microscope (AFM) to 
measure the Casimir force between a small sphere and a fiat plate (140, 163, 164, 
165). Their most recent measurement used a 191-//m-diameter polystyrene sphere 
that was glued to a 320-/xm-long AFM cantilever. The cantilever plus sphere and 
a 1-cm-diameter optically polished sapphire disk were coated with 87 nm of gold, 
with a measured surface roughness of 1.0 ± 0.1 nm. The disk was placed on a 
piezoelectric tube with the sphere mounted above it, as shown in Figure 10. The 
cantilever flex was measured by reflecting laser light from the cantilever onto split 
photodiodes. The force scale was calibrated electrostatically by applying a zb3-V 
potential difference between the sphere and disk at a separation of 3 //m. The 
force difference between the -|-3 V and —3 V applied potentials was used to de- 
termine the residual potential difference between the disk and sphere when their 
external leads were grounded together: 3 ± 3 mV. The force between the sphere 
and disk was measured for separations ranging from 400 nm to contact. It was 
found that the surfaces touched when their average separation was 32.7 ±0.8 nm. 
This was attributed to gold crystals protruding from the surfaces. The measured 
forces were compared to the expected Casimir force for separations of 62-350 nm 
and agreement to within 1% was found (165). 
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The record for measuring the Casimir force at the closest separation is held by 

Ederth (166) at the Royal Institute of Technology in Stockholm, who measured 

the force between crossed cylindrical silica disks with diameters of 20 mm. A 

template-stripping method (167) was used to glue 200-nm layers of gold, with an 

rms surface roughness of < 0.4 nm, to the silica disks. The gold surfaces were then 

coated with a 2.1-nm-thick layer of hydrocarbon chains to prevent the adsorption 

of surface contaminants and the cold-welding of the gold surfaces upon contact. 

One cylindrical surface was attached to a piezoelectric stack and the other to a 

piezoelectric bimorph deflection sensor that acted as a cantilever spring. The two 

surfaces were moved toward one another starting at a separation > 1 //m, where 

the Casimir force was less than the resolution of the force sensor, and ending 

at a separation of 20 nm, at which point the gradient of the Casimir force was 

comparable to the stiffness of the bimorph spring, causing the surfaces to jump 

into contact. The stiffness of the bimorph sensor was calibrated by continuing 

to move the piezotube another 200-300 nm while the surfaces were in contact. 

The absolute separation between the surfaces was found by fitting the measured 

force curve to the expected Casimir signal (plus electrostatic background, which 

was found to be negligible) with the absolute separation as a fit parameter. At 

contact, the surfaces compressed by w 10 nm. The measured force was compared 

to the expected Casimir force over the range of separations from 20 to 100 nm 

and an agreement to better than 1% was found. 

4.4.2 Signal-to-noise and background considerations 

The signal-to-noise ratio for Casimir-force measurements as tests of the ISL may 
be improved by using more sensitive force probes, using thicker metallic coat- 
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ings on the test bodies, and operating at lower temperatures. Nonetheless, the 

dominant limitation for interpreting the measurements as tests of the ISL is un- 
derstanding the Casimir-force background to high accuracy. There is a growing 
literature on the corrections that must be applied to the Casimir force calculated 
for smooth, perfect conductors at zero temperature (Equations 35 and 36). The 
dominant corrections are for finite temperature, finite conductivity, and surface 
roughness. Corrections for finite temperature are important for test-body sepa- 
rations d > 1 iJ,m. For the Lamoreaux experiment, the finite-temperature correc- 
tions at l-/xm and 6-/xm separations were 2.7% and 174% of the zero-temperature 
Casimir force, respectively (168). A number of authors have considered the ef- 
fects of finite conductivity on the temperature correction (169, 170, 171), and 
results believed to be accurate to better than 1% were obtained. The correction 
to the Casimir force for the finite conductivity of the metallic surfaces is of or- 
der 10% at d = 1 /xm and grows with smaller separations. Finite-conductivity 
corrections using a plasma model for the dielectric function of the metal give 
the correction as a power series in Xp/d, where Ap is the plasma wavelength 
of the metal (172, 173, 174). Corrections have also been obtained using optical 
data for the complex dielectric function (175, 165, 176, 177). Surface roughness 
of the test bodies contributes a correction to the Casimir force that can be ex- 
pressed as a power series in h/d, where ^ is a characteristic amplitude of the 
surface distortion (178, 179, 180, 174). For stochastic distortions, the leading- 
order surface-roughness correction is 6{h/d)'^, which is less than 1% of the Casimir 
force at closest separation in the experiments of Ederth and the Riverside group. 
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4.4.3 Results 

Constraints on Yukawa interactions with ranges between 1 nm and 10//m, shown 
in Figure 9, have been extracted from the Casimir-force measurements of Lam- 
oreaux (168, 151), Ederth (181), and the Riverside group (182, 183, 7, 184). 
Figure 9 also shows constraints at even smaller ranges obtained from earlier van 
der Waals-force experiments (185). It should be noted that most of these con- 
straints were obtained by assuming that a Yukawa force could not exceed the 
difference between the measured force and the predicted Casimir effect. To be 
rigorous, the raw data should be fitted simultaneously with both Casimir and 
Yukawa forces, which should lead to significantly less stringent limits on |q:|. 
Deviations from Newtonian gravity in this region that follow a power law (Equa- 
tion 3) are constrained more strongly by the much more sensitive longer-range 
gravity experiments discussed above (7). 

4-5 Astronomical Tests 

A summary of constraints on Yukawa interactions with A > 1 mm may be found 
in Figure 2.13 of the 1999 review by Fischbach &; Talmadge (14), which we 
reproduce in part in our Figure 4. Since the publication of Reference (14), the 
constraints for A < 1 cm have been substantially improved, as discussed above. 
The constraints at larger ranges from laboratory, geophysical, and astronomical 
data (see Figure 4) are essentially unchanged. The astronomical tests provide the 
tightest constraints on a. These are typically based on Keplerian tests comparing 
G{r)MQ values deduced for different planets. However, the tightest constraint 
comes from lunar-laser-ranging (LLR) studies of the lunar orbit. Because this 
result may improve significantly in the next few years, we give some details of 
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the measurement here. 

The LLR data consist of range measurements from telescopes on Earth to 
retrorefiectors placed on the Moon by US astronauts and an unmanned Soviet 
lander. The measurements, which began in 1969, now have individual raw range 
precisions of about 2 cm and are obtained from single photon returns, one of which 
is detected for roughly every 100 launched laser pulses (186). The vast majority 
of the data come from sites in Texas (187) and in southern Prance (188). The 
launched laser pulses have full widths at half maximum of about 100 ps; the return 
pulses are broadened to about 400 ps because the reflector arrays typically do 
not point straight back to Earth owing to lunar librations. The launch-telescope- 
to-lunar-retroreflector ranges have to be corrected for atmospheric delay, which 
is computed from the local barometric pressure, temperature, and humidity. For 
the Moon straight overhead, the range correction at the Texas site is about 2 m. 
The dominant uncertainties in converting raw range measurements into separa- 
tions between the centers of mass of the Earth and the Moon come from tidal 
distortions of the Earth and Moon and atmospheric and ocean loading of the 
Earth. The current model, using the entire world data set, gives an uncertainty 
of about 0.4 cm in the important orbit parameters. 

The most sensitive observable for testing the ISL is the anomalous precession 
of the lunar orbit. If the Moon were subject only to a central Newtonian 1/r po- 
tential from the Earth, the lunar orbit would not precess. The orbit does precess 
due to the Earth's quadrupole field and perturbations from other solar-system 
bodies, as well as from the small general relativistic geodetic precession and pos- 
sibly also from a Yukawa interaction; the conventional sources of precession must 
be accounted for to obtain the anomalous Yukawa precession rate. Ignoring terms 
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of order e^, where the Moon's eccentricity is £ = 0.0549, the anomalous Yukawa 
precession rate Su is (14) 



where a; = 27r radians/month and a is the mean radius of the Moon's orbit. The 
constraint on a (A) is tightest for A = a/2 and falls off relatively steeply on either 
side of A = a/2. The current LLR 2a upper limit on Su is 270 /xas/y; this follows 
because the observed precession of about 19.2 mas/y agrees with the general 
relativistic prediction to (— 0.26ib0.70)%, where the error is "realistic" rather than 
"formal" (the error quoted in Reference (189) should be doubled; J. Williams, 
private communication 2003). We conclude that at 95% CL, Slo/lo < 1.6 x 10~^^; 
the corresponding LLR constraint is shown in Figure 4. 



5. 1 Summary of Experimental Results 

Because gravity is intimately connected to the geometry of spacetime, ISL tests 
could provide very direct evidence for the existence of extra space dimensions. 
In addition, ISL tests are sensitive to the exchange of proposed new low-mass 
bosons. A variety of theoretical considerations hint that new effects may oc- 
cur at length scales between 10 jim. and 1 mm. This circumstance, as well as 
the urge to explore unmapped territory, has motivated the development of new 
experimental techniques that have produced substantial improvements in con- 
straints on theories. The overall slope of the experimental constraints shown in 
Figures 5, 9, and 4 reflects the rapidly decreasing signal strength of a new inter- 
action as its range decreases. At gravitational strength (a = 1 in Figure 5), the 
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ISL has been verified down to a distance A = 200 jim. At length scales between 

20 nm and 4 mm, many square decades in Yukawa-parameter space have been 

ruled out. These results have eliminated some specific theoretical scenarios, but 

many other interesting ideas are still viable because their predicted eff'ects lie 

somewhat below the current experimental limits. 

5. 2 Prospects for Improvements 
5.2.1 Short-range tests of the ISL 

To make a gravitational-strength {a = 1) ISL test at a 20-/xm length scale re- 
quires an increase in the background-free sensitivity of at least a factor of 10^. 
Fortunately, such an increase is possible, although it will require years of devel- 
opment. 

The Eot-Wash group are currently running a new apparatus that features a 
pendulum/attractor system having 22-fold rotational symmetry with 44 thinner, 
smaller-diameter holes. The pendulum ring and attractor disk are made from 
denser materials (copper and molybdenum, respectively). Noise has been im- 
proved by a factor of six. The closest attainable separation has been reduced 
by a factor of two by adding a passive "bounce" -mode damper to the fiber- 
suspension system, and the thickness of the electrostatic shield has been reduced 
to 10 nm. Figure 3 shows the spectral density of the torque signal from this 
apparatus. This instrument should probe Yukawa forces with \a\ = 1 for ranges 
down to A = 60 fxm. In principle, it is possible to use a low-frequency torsion 
balance in a different mode, one that measures the attraction between two flat 
plates (J.G. Gundlach, private communication). This would provide a null test 
with a sensitivity that scaled as }?e~^/^ rather than as \^e~^/^ in the partial-null 
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experiments. 

The Colorado group plans to optimize their geometry and to use a Washington- 
style electrostatic shield to attain closer separations. This could improve their 
limits between 10 fj,m and 50 /xm by at least an order of magnitude. In the long 
run, both groups could run at liquid-helium temperatures, which will give lower 
noise, not only from the decreased ksT factor, but also from the expected in- 
crease in the Q factor of the torsion oscillator. Newman (190) found that the Q 
factor of a torsion fiber has two components. One is temperature-independent but 
amplitude-dependent (this is already negligible in the Eot-Wash instrument be- 
cause of the small amplitudes employed) and the other is temperature-dependent 
and amplitude-independent. 

The microcantilever application exploited by the Stanford group has not yet 
attained its full potential. Presumably, lessons learned in this pioneering exper- 
iment will reduce the backgrounds and allow the experimenters to exploit their 
inherent sensitivity to new very small forces. Because corrections to the ideal- 
ized Casimir force can be large and depend on properties of the test bodies that 
are troublesome to quantify, it may be difficult to compare Casimir-force exper- 
iments to theory at an accuracy much better than 1%. The finite-conductivity 
corrections depend on the dielectric properties of the actual metallic coating of 
the test bodies, which may differ somewhat from bulk dielectric properties used 
in the calculation. As the experimental precision improves, parameters associ- 
ated with the conductivity correction (such as Ap) may need to be included as 
adjustable parameters in fitting the measured force-versus-distance curves. The 
surface-roughness correction should consider distortions over length scales larger 
than are easily accessible by AFM scans, and it may be necessary to vary the 
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roughness parameters as well. Both corrections scale as inverse powers of the 

separation, d, as do the corrections for residual electric potential and patch ef- 
fects. Compounding the problem of multiple corrections with similar distance 
dependences is the uncertainty in the absolute separation of the test bodies. The 
Casimir force depends on do + dr rather than on d, where dr is the relative 
displacement of the test bodies between force measurements (which can be accu- 
rately measured) and do is the absolute separation at the origin of the relative 
scale (which is difficult to determine accurately) . Including do as a fit parameter 
allows other short-distance parameters to vary (166) without affecting the fit at 
large distances, where the fractional error on the force measurements is larger. 
It is unlikely that the next few years will bring large improvements in Yukawa 
constraints from Casimir-force experiments. 

5.2.2 Long-range tests of the ISL 

Because any change in orientation of the Moon's ellipticity grows linearly with 
time, even with data of constant precision the LLR constraint should improve in 
proportion as the data span increases (assuming that the modeling of conventional 
precession sources is not a limiting factor). New LLR projects should improve 
the raw range precision by an order of magnitude, bringing the precision into the 
range needed to test the "high-pass" gravity model (6). For example, APOLLO 
(191) will exploit a 3.5-m telescope at an elevation of 2780 m and subarcsecond 
image quality. This instrument should receive several returned photons per laser 
shot, giving a data rate about 10^ times greater than existing facilities. It is 
expected that more precise data will lead to corresponding improvements in the 
modeling. 
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Ranging to other planets is necessary in order to probe longer length scales 

effectively. This is currently done using radar (which is limited by the absence of 
a well-defined "target" on the planet) or else microwave signals transmitted by 
orbiting spacecraft (which are limited by uncertainties and the finite timespan 
of the orbits). Furthermore, the accuracy of microwave ranges is limited by 
propagation delay in the interplanetary solar plasma. It is impractical to laser- 
range to passive reflectors on other planets (if they could be placed) because 
the returned signal falls as However, recent developments in active laser 

transponders, whose sensitivity falls as make it practical to place such a 

device on Mars and ultimately achieve range precisions of a few centimeters (192). 
This would yield several interesting new gravitational measurements, including 
an improved test of the strong equivalence principle (193), which provides one 
of the best limits on massless gravitational scalar fields, as well as tests of the 
ISL that would give interesting constraints on the quasilocalized gravity model 
of Reference (128). 

ISL tests at scales larger than the solar system typically rely on uncertain 
astrophysical models. But Will (5) notes that the proposed LISA space-based 
interferometer could test a pure Yukawa potential at a scale of 5 x 10^^ m by 
studying distortions of the gravitational waveform from an inspiraling pair of 
1O^M0 compact objects. 

5.3 What if a Violation of the Law Were Observed? 

Suppose that future experiments revealed a violation of the ISL at short length 
scales. Of course, one would try to tighten the constraints on its range and 
strength by performing tests using instruments with varying length scales. But 
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a new question immediately arises: Is the new physics a geometrical effect of 

extra dimensions or evidence for exchange of a new boson? This can be decided 
by testing whether the short-range interaction violates the equivalence principle: 
Boson exchange generically does not couple to matter in a universal manner and 
therefore appears as a "violation" of the equivalence principle, whereas geomet- 
rical effects must respect the principle. Kaplan & Wise (102) estimated that the 
equivalence-principle- "violating" effect from dilaton exchange is 0.3%. 
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Figure 1: Schematic diagram of the 10-hole pendulums and rotating attractors 
used in the two experiments of Hoyle et al. (71, 148, 149). The active components 
are shaded. 
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Figure 2: Top; Torques measured in the first experiment of Hoyle et al. as a 
function of pendulum/attractor separation. Open circles are data taken with the 
lower attractor disk removed and show the effect of uncancelled gravity. Smooth 
curves show the Newtonian fit. Bottom: Residuals for the Newtonian fit. The 
solid curve shows the expected residual for a Yukawa force with a = 3 and 
A = 250 /xm. 
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Figure 3: Spectral density of the torque signal in the 22-fold symmetric exper- 
iment of the Eot-Wash group. The peaks at 8.5 and 17 to are gravitational 
calibrations; the fundamental and first three overtones of the short-range signal 
are at 22, 44, 66, and 88 to. The smooth curve shows the thermal noise computed 
using Equation 33. 
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Figure 4: 95%-confideiice-level constraints on ISL-violating Yukawa interactions 
with A > 1 cm. The LLR constraint is based on the anomalous perigee precession; 
the remaining constraints are based on Keplerian tests. This plot is based on 
Figure 2.13 of Reference (14) and updated to include recent LLR results. 
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Figure 5: 95%-confidence-level constraints on ISL-violating Yukawa interactions 
with 1 < A < 1 cm. The heavy curves give experimental upper Hmits (the 
Lamoreaux constraint was computed in Reference (151)). Theoretical expecta- 
tions for extra dimensions (56), moduli (101), dilaton (102), and radion (83) are 
shown as well. 



78 



ADELBERGER, HECKEL k NELSON 




Tuning Block 



Figure 6: Schematic diagram of the instrument used by Long et al. (141). 
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Figure 7: Data from the experiment of Long et al. (141) showing the two quadra- 
ture signals from the torsion oscillator. 




Figure 8: Schematic diagram of the instrument used by Chiaverini et al. (152). 
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Figure 9: Constraints on ISL-violating Yukawa interactions with Inm < A < 1/im 
adapted from Reference (7). As discussed in the text, these upper hmits, ex- 
tracted from Casimir-force measurements, are not as rigorous as those in Fig- 
ures 5 and 4. 
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Figure 10: Schematic diagram of the Casimir-force apparatus used in Refer- 
ence (165). 



